Constraints on neutrino — dark matter interactions from cosmic microwave 
background and large scale structure data 
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We update a previous investigation of cosmological effects of a non-standard interaction between 
neutrinos and dark matter. Parameterizing the elastic-scattering cross section between the two 
species as a function of the temperature of the universe, the resulting neutrino-dark matter fluid has 
a non-zero pressure, which determines diffusion-damped oscillations in the matter power spectrum 
similar to the acoustic oscillations generated by the photon-baryon fluid. Using cosmic microwave 
background data in combination with large scale structure experiment results, we then put con- 
straints on the fraction of the interacting dark matter component as well as on the corresponding 
opacity. 
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I. INTRODUCTION 



Both the existence and the energy content of dark 
matter in the universe have been firmly established 
by several observations (see [lj for a review), yet, after 
decades of research its nature is still unknown. The 
favored candidate for dark matter has three main fea- 
ture namely, it seems to be a long-lived, cold and col- 
lisionless species. However, it is possible to relax some 
of these assumptions and in fact, several different dark 
matter candidates have been proposed through years 
( @j 9 ) an d it is not even clear if only one species is 
responsible for the whole dark matter in the universe. 
An interacting component has been proposed by many 
authors; couplings with a light scalar or pseudoscalar 
boson, as in the Majoron model has been investigated 
in 0, H, H, and interaction with neutrinos and more 
generally with the electromagnetic plas ma have been 
considered for example, in [cTl9l. [ToL llll. [l2|. 
In this paper we study the scenario where at least a 
fraction of dark matter is coupled to neutrinos. We 
parameterize the neutrino-dark matter interaction in 
terms of the opacity P = {o~d m ~u\v\) /rridm, the ratio of 
the thermal averaged dark-matter-neutrino scattering 
cross section to the mass of the dark matter particle, 
where P = Q a~ 2 with Q a constant and a the scale 
factor. We will briefly review the underlying theoret- 
ical basis of this choice in the next Section. 
The main observable consequence of such a coupling 
is a pattern of oscillations in the matter power spec- 
trum, which can in principle mimic the baryon acous- 
tic oscillations observed in the large scale clustering 
of galaxies [H], [HI- The cosmic microwave back- 
ground (CMB) radiation is also affected by this in- 
teraction. Indeed, for sufficiently small scales, the 
temporal growth of a dark matter perturbation is low- 
ered because of the coupling with neutrinos, and this 
determines a small enhancement of the small scale 
peaks at a level which depends on both the strength of 



the interaction and the percentage of interacting dark 
matter. This effect can be also described in terms of 
a lowering of the neutrino "viscosity parameter" c 2 is 
[l|| from its standard value c 2 vis = 1/3. 

This paper updates and improves the results of a 
previous work [131 ]. In fact, we present the results of a 
likelihood analysis based on a Monte Carlo Markov 
Chain sampling of the full cosmological parameter 
space, which allows us to take into account for possi- 
ble degeneracies with the new parameters introduced 
in the analysis. To this aim, we use the latest CMB 
datasets available from the WMAP experiment [ijj 
in combination with large scale structure data from 
SDSS. 

The paper is organized as follows. In Section II we 
briefly discuss possible theoretical models for the dark 
matter-neutrino interaction and corresponding opac- 
ity due to scatterings. We refer to [13[ for a more 
detailed discussion of these topics. In Section III we 
explain the method considered in the analysis and the 
data used to constrain both the strength of the in- 
teraction and the fraction of interacting dark matter. 
Finally we give our results and conclusions. 



II. NEUTRINO-DARK-MATTER 
SCATTERINGS 

Neutrino - dark matter interactions can be mod- 
elled in several ways, depending on the spin of the 
dark matter candidates as well the specific coup ling 
lagrangian. This issue is covered in details in LL3| 
to which we refer for more details, while we summarize 
in this Section the main general features of such inter- 
actions and report the thermal averaged cross section 
times velocity as function of the scale factor. The lat- 
ter, denoted by a is normalized to unity today, and 
scales as the inverse neutrino temperature. 

As was shown in [l3j , if scatterings are supposed to 
affect the structure formation on large scales, it is un- 
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been also considered in the literature, where dark mat- 
ter relic density today is the remnant of an initial par- 
ticle - antiparticle asymmetry produced in the early 
universe, thus similar to the mechanism of baryogene- 
sis which leads to a baryon density today much larger 
than what is expected by freezing of strong interac- 
tions alone [HI, This possibility may also ac- 
count for the intriguing similarity between the values 
of Qi, and Qdm we observe today, which differ by a fac- 
tor five only, a feature which may call for a common 
mechanism for their formation. In the following we 
therefore, consider the case of a non self-conjugated 
fermion or scalar particle ip with a conserved global 
U(l) charge, at least in the low energy scale regime, 
say below their mass scale, which also corresponds to 
the relevant stages for structure formation we are in- 
terested in. 

In case of a spin zero species, the interaction la- 
grangian can be chosen as a Yukawa term 



hFRVLijj + h.c. 



(1) 



with F a spinor field, or via coupling with an inter- 
mediate vector-boson field 

Ant = z 5 (W-W)^+fl 2 f^^" 
+ guVL'fvi.Uft. (2) 
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FIG. 1: A perturbation of interacting dark matter (red 
line) of wavenumber k = 0.81 ZiMpc" 1 is plotted against a 
perturbation of non-interacting dark matter (black dotted 
line) in the upper panel. The effects of this interaction 
are clearly seen in the angular power spectum of the CMB 
(middle panel) and on the matter power spectum (lower 
panel) . We have chosen a very large value of Q to magnify 
the effect on the main cosmological observables. 



In both cases F and U fields will be assumed to have 
mass larger than ip, to prevent fast ip decay at tree 
level, see [l3j |. 

In case of spin 1/2 dark matter, one has 



Ant = hip R v L F + h.c. 
with F a scalar field or finally, 



(3) 



Ant = gicLtp^IpL + CRlpj^lp^Un 

+ guVLl^VLU^. (4) 

We notice that, as we will see in the next Section, 
bounds on dark matter - neutrino scattering cross sec- 
tion from CMB and large scale data, correspond to a 
mass scale for ip as well as the other field involved in 
the interaction lagrangian (F or U) larger than MeV. 
Thus, we will assume this lower bound in the follow- 
ing, m dm > MeV, rriF.u > MeV. 

It is quite easy to compute the thermal averaged 
scattering cross section corresponding to these inter- 
action terms in the non-relativistic limit for the dark 
matter particle ip and F or U . For example in the 
case of Eq. ^ or © and (f2]) one gets, respectively 



likely that the relic abundance of these non standard 
dark matter candidates formed via the well known 
freezing phenomenon, as the pair annihilations would 
be quite efficient till recent times in depleting their 
number density, unless they represent a too light 
species. On the other hand, a different scenario has 
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The first result holds for mp ^ mdm, otherwise cross 
section takes a constant value, as for the well known 
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TABLE I: Mean values and marginalized confidence 
levels for the cosmological parameters of interest for 
WMAP+SDSS. When not otherwise stated, all constraints 
are at 68% c.L. 



Parameter 


WMAP+SDSS 
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Thomson cross section. This special case has been 
considered in details in (TH, and will not be worked 
out here, as it is quite fine-tuned. Similar results hold 
for a spinorial dark matter interacting via an interme- 
diate gauge boson U 
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We see that for both Yukawa-type or intermediate 
gauge boson interaction, the opacity assumes a simple 
expression, of the form 



1 



M 4 m dl 



rp2 



1 V 



(0)2 



= Q a 



(8) 



with M a mass scale and Q a constant parameterizing 
the strength of interaction 
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cm 2 MeV _1 , (9) 



where we have used the value of neutrino temperature 
today Tj 0) . 



III. ANALYSIS 

In the early universe, the opacity P is strong 
enough to determine acoustic oscillations in the cou- 
pled neutrino-dark matter fluid, just like the usual 
electromagnetic force does with the photon-baryon 
fluid before recombination. In order to include a 
neutrino-dark matter interaction, we modify the stan- 
dard Eulero equation for dark matter, resulting in an 
analogue of that for baryon matter 
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FIG. 2: Constraints in the cross section parameter Q 
versus Q,xh 2 plane from the combination of WMAP and 
SDSS data. Large values for the neutrino-dark matter in- 
teraction are always compatible with the standard picure, 
provided that the energy content of the interacting dark 
matter is small enough. 



where momentum conservation in scattering processes 
has been accounted for. With 9 we denote the velocity 
perturbations and the product an dm <?dm-v quantifies 
the scattering rate of neutrinos by dark matter. As in 
the case of the baryon-photon interaction, we can ne- 
glect both the shear contribution k 2 a and the term 



-"dm 



k 5, where c dm is the sound speed of the dark- 
matter fluid. The hierarchy of equations for neutrinos 
has been modified accordingly. 

As we can see from the upper panel of Fig. 1 , any inter- 
acting dark matter perturbation of sufficiently small 
scale will be affected by the coupling with neutrinos. 
When it enters the horizon, the mode will experience a 
series of acoustic oscillations while coupled with neu- 
trinos, rather than growing sligthly in radiation dom- 
inated epoch and then more rapidly in matter domi- 
nated epoch. This determines a pattern of oscillations 
in the matter power spectrum (lower panel) as well as 
an enhancement of the peaks in the angular power 
spectrum (middle panel). 

We make use of the publicly available Markov Chain 
Monte Carlo (MCMC) package CosmoMC [1| with 
a convergence diag nostic based on the Gelman and 
Rubin statistic [2l| to put constraints on the cosmo- 
logical parameters. We consider a fiat ACDM model 
described by the following set of nine cosmological pa- 
rameters: 

{(j b ,Lj c ,u)x,® s ,T,n s ,\og[10 w A s },Q,Asz} , (H) 

where u>b (= Qbh 2 ), uj c (= Q c h 2 ), (= ^xh 2 ) are 
the baryon, standard collisionlcss cold dark matter 
and interacting cold dark matter densities relative to 
the critical density respectively, 9 S is the ratio of the 
sound horizon to the angular diameter distance at de- 
coupling, t the optical depth to reionization and A s 
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and n s are the amplitude of the primordial spectrum 
and the spectral index. They are both evaluated at 
k = 0.002 Mpc -1 . The parameter Q (measured in 
units of cm 2 MeV -1 ) quantifies the strength of the 
neutrino-dark matter interaction, as discussed in the 
previous section. Finally, we take into account for 
the contribution to the power spectrum by Sunyaev- 
Zcl'dovich fluctuations 1231 by adding the predicted 
template spectrum from |23l |. multiplied by the con- 
stant amplitude Agz- We use data and likelihood 
code from the WMAP team's 5 year release [13, HI 
(both temperature and polarization) in combination 
with the SDSS LRG data release [25j . 
In Fig. 2 we plot the combined constraints on both 
the parameter Q and the energy content in interacting 
dark matter. As expected, there is a clear correlation 
between these two parameters as the lower the energy 
content in interacting dark matter the weaker the con- 
straint on the strength of the interaction. It has been 
shown in (l3j that the baryon acoustic oscillations due 
to the coupled photon-baryon plasma and seen in the 
large scale clustering of matter [l4|, EH can mimic the 
coupled neutrino-dark matter fluid in a universe with 
a lower energy content in baryon. seriously affecting 
constraints on fif,/i 2 derived by large scale structures 
data only. We find here that when a full parame- 
ter space analysis is performed using additionaly data 
from the CMB, we always recover the standard con- 
straints for the baryon density. Of course the energy 
content in standard dark matter is smaller than the 



usual amount simply because we are splitting all the 
dark matter in an interacting and a non interacting 
component. All the other values for the cosmological 
parameters are compatible with the standard values 
at la. 



IV. CONCLUSIONS 

In this paper we have investigated the effects of a 
non-standard interaction between neutrinos and dark 
matter using the to date available data from WMAP 
and SDSS Collaborations on cosmic microwave back- 
ground and large scale structures. These effects are 
the presence of diffusion - damped oscillations in the 
matter power spectrum and a small enhancement of 
the peaks in the cosmic microwave background angu- 
lar power spectrum. Parameterizing the opacity pa- 
rameter as Q a~ 2 with Q a constant, which holds in 
a large class of models for neutrino - dark matter in- 
teractions, we found that the bound on Q becomes 
stronger as the percentage Qxh 2 of interacting dark 
matter component grows. Typical upper bound on Q 
is of the order of 10" 42 ™ 10" 40 cm 2 MeV" 1 , while the 
bound becomes weaker for very low values of fix] we 
note that cosmological constraints on neutrino masses 
can be affected because of possible degeneracies. 

This research was funded by NSF CAREER AST- 
0605427 and by LANL IGPP-08-505. 
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